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Sorin DragomiiQ and Stefano Montaldd 

Abstract. We study subelliptic biharmonic maps i.e. smooth 
maps <p : M — > N from a compact strictly pseudoconvex CR man- 
ifold M into a Riemannian manifold N which are critical points 
of the energy functional E 2 , b (<P) = \ f M \\r b (<f>)\\ 2 6 A {d6) n . We 
show that <j) : M —¥ N is a subelliptic biharmonic map if and only 
if its vertical lift <f> f : C(M) -4 JV to the (total space of the) 
canonical circle bundle S 1 — > C(M) — ^» M is a biharmonic map 
with respect to the Fefferman metric Fg on C(M). 

1. Introduction 

Biharmonic maps were introduced by J. Eells & L. Lemaire, [18], as 
critical points <fi G C°°(M, N) of the bienergy functional 



£ 2 (0) = 5 / 



|r(0)|| 2 rf^ 



where M and are Riemannian manifolds and r(<p) is the tension field 
of : M — > N. Biharmonic maps were further investigated by G. Jiang, 
[22] . who derived the first and second variation formulas for E 2 (<fi) and 
gave several applications to the geometry of the second fundamental 
form of a submanifold in a Riemannian manifold. In the last decade 
a large amount of work has been devoted to biharmonic maps with 
particular attention to the constructions and classifications of proper 
biharmonic maps and proper biharmonic submanifolds, see e.g. [2]- [3], 

A program aiming to extending results on nonlinear elliptic systems 
of variational origin to the hypoelliptic case was started by J. Jost & 
C-J. Xu, [23J. Given a Hormander system of vector fields {Xi, • • • , X p } 
on an open set U C M m and a map <p G C°°(U, N) the function 

1 P 

e(0) = -^X a (0^)X a (0 fc )(^o0) 
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is globally defined and generalizes the ordinary energy density of 0. 
Here hjk is the Riemannian metric on N in a local coordinate system 
(V, yi) and <p> = yi o 0. A subelliptic harmonic map is a critical point 
G C°°(U, N) of the functional (cf. [23]) 

Ex{4>) = [ e(0) dx 
Jn 

where Q C M m is a bounded domain such that Q C U. The Euler- 
Lagrange equations of the variational principle S Ex{4>) = are 

p 

-h^ + (r jk o 0) ^x a (0^)x a (0 fc ) = o 

a=l 

where if = Y7 a =i -^a-^-a is the Hormander operator and T l - k are the 
Christoffel symbols of hjk- Subelliptic harmonic maps were recognized 
(cf. E. Barletta & S. Dragomir & H. Urakawa, [8]) as the local man- 
ifestation of pseudoharmonic maps i.e. critical points G C°°(M, N) 
of the functional 

ExM = lf trace Ge (U H <j>*h) 9 A (d6) n . 

Here M is a compact strictly pseudoconvex CR manifold, of CR dimen- 
sion n, and 6 is a contact form on M such that the Levi form Gq is pos- 
itive definite. Also h is the Riemannian metric on iV and n# <p*h is the 
restriction of the bilinear form (p*h to the Levi, or maximally complex, 
distribution H(M). The Euler-Lagrange equations of 5 Ei^((j)) = 
may be written as Tb(0) = where the field Tfc(0) G C oc (0~ 1 T(A^)) is 
locally given by 

2n 

n(<py = A b 0* + ^x a (0^)x a (0 fc ) (r jk o 0) . 

a=l 

Here {X a : 1 < a < 2n} is a local GVorthonormal frame of H(M). 
Also Afe is the sublaplacian i.e. the second order differential operator 
given by 

A b u = div (V H u) , u G C°°(M). 

The divergence operator is meant with respect to the volume form 
9 A (d6) n while the horizontal gradient is given by V H u = n^Vw and 
geCVu,X) = X(u) for any X G 3£(M). CR manifolds occur mainly 
as boundaries M = dfl of domains Q C C n+1 and boundary values 
of Bergman-harmonic maps $ : Q — > N may be shown (cf. [16]) to 
be pseudoharmonic provided $ has vanishing normal derivatives (thus 
motivating our use of the index b for "boundary" analogs to geometric 
objects such as the tension field t($), the second fundamental form 
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etc.). The similar boundary behavior of biharmonic maps from 
a strictly pseudoconvex domain Q C C n+1 endowed with the Bergman 
metric is unknown. 

The approaches in [23] and |8J overlap partially, as follows. For any 
local Gg-orthonormal frame {X a : 1 < a < 2n} in H(M) defined 
on a local coordinate neighborhood ip : U — >■ R 2n+1 the push-forward 
{(dip)X a : 1 < a < 2n} is a Hormander system on (p(U) and A& = —H. 
However in [23] and in general in the theory of Hormander vector fields 
the Euclidean dimension m is arbitrary (as opposed to m = 2n + 1 in 
the CR case) , the vector fields forming the given system are allowed to 
be linearly dependent (at particular points), and the formal adjoints 
X* are meant with respect to the Euclidean metric on Q (rather than 
the Webster metric - a non flat Riemannian metric springing from the 
given CR structure in the presence of a contact form). It is therefore a 
natural problem, within J.Jost & C-J. Xu's program, to study critical 
points (f) G C°°(M, N) of the functional 



These are referred to as subelliptic biharmonic maps. We give a geomet- 
ric interpretation of subelliptic biharmonic maps in terms of ordinary 
biharmonic maps from a Lorentzian manifold (the total space of the 
canonical circle bundle S 1 — > C(M) — > M endowed with the Feffer- 
man metric Fg, [1 TJ ) . The paper is organized as follows. In §2 we 
consider ordinary biharmonic maps from a Fefferman space-time. The 
rough Laplacian (a degenerate elliptic operator appearing in the prin- 
cipal part of the subelliptic biharmonic map system) is discussed in § 3. 
The first variation formula for the functional E 2i b '■ C°°(M, N) — > R is 
derived in §4 while the main result (cf. Theorem [1] below) is proved in 
§ 5. A few open problems are outlined in § 6. 

2. Biharmonic maps from Fefferman space-times 

Let (M, Ti o(M)) be a compact orientable CR manifold, of CR di- 
mension n, where T\$(M) is its CR structure. Let us assume that 
M is strictly pseudoconvex and 6 a contact form on M such that the 
Levi form Gg(X, Y) = (dO)(X,JY), X, Y e H(M), is positive defi- 
nite. Here H(M) = Re{T lj0 (M) ©T ,i(M)} is the Levi distribution 
of M and J : H(M) H(M), J(Z + Z) = i(Z -Z), Z e T lfi (M), 
its complex structure (z = \f—l). For all needed notions of CR and 
pseudohermitian geometry we rely on [17] . Let T be the Reeb vector of 
(M, 9) i.e. the nowhere zero globally defined tangent vector field field 
transverse to H(M) determined by 9(T) = 1 and TJ d9 = 0. Let V be 
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the Tanaka- Webster connection of (M,9) (cf. Theorem 3.1 in [T7j, p. 
25, for the axiomatic description of V). 

Let S 1 — > C(M) — > M be the canonical circle bundle over M and 
Fg the Fefferman metric on C(M) (cf. J.M. Lee, [23], or Definition 2.15 
in [IT], p. 128) i.e. the Lorentzian metric on C(M) given by (cf. (2.30) 
in [HI, p. 127) 

(1) Fg = n*Gg + 2 (n*9) a, 

(2) "~ {« + «* (*< - \ - ^ O) } ■ 

As to the notation in (JTJ)- (J2J) , given a local frame {T a : 1 < a < n} 
of Tifi(M) we set g a g = Gg(T a ,T-p) (with T„ = T a ). Also are 
the connection 1-forms of V with respect to {T a : 1 < a < n} i.e. 
VT/3 = Up g) T a . Moreover p = g a/3 R a p is the pseudohermitian scalar 

curvature of (M, 9) (cf. e.g. [TF], p. 50). The (0,2)-tensor field G e 
in ([I]) is obtained by extending the Levi form Gg to a degenerate form 
defined on the whole of T(M). By definition Gg = Gg on H(M)®H(M) 
and Gg(X, T) = for any X G T(M). By a result of C.R. Graham, [20], 
the (globally defined) 1-form a is a connection form in the canonical 
circle bundle. Let X t G X(C(M)) be the horizontal lift of X G X(M) 
with respect to a. If S G 3£(C(M)) is the tangent to the S 1 action 
then — S is a nowhere vanishing globally defined timelike vector 
field hence (C(M), Fg) is a time-oriented Lorentzian manifold, referred 
to as the Fefferman space-time. Fg was discovered by C. Fefferman, 
[T5] , in connection with the study of the boundary behavior of the 
Bergman kernel of a bounded domain flcC" (as a Lorentzian metric 
on C(dQ) ~ dQ x S 1 ). An array of popular nonlinear problems arise 
from Fg e.g. the CR Yamabe problem (cf. |17] . p. 159-160) is the 
projection via 7r : C(M) — > M of the ordinary Yamabe problem for Fg. 
While the principal part of the Yamabe equation on C (M) is the wave 
operator □ (hence the Yamabe equation on C(M) is not elliptic) the 
principal part of the projected equation is the sublaplacian A b (hence 
subelliptic theory applies, cf. [IT], p. 176-210). Also pseudoharmonic 
maps from strictly pseudoconvex CR manifolds may be characterized 
as base maps of 5' 1 -invariant harmonic maps from (C(M),Fg), thus 
suggesting that Theorem [1] below should hold. 

Several basic facts in harmonic map theory are known to extend in a 
straightforward manner from the Riemannian to the semi-Riemannian 
setting (cf. [1], p. 427-452) e.g. a smooth map of semi-Riemannian 
manifolds has a well defined tension tensor field. In particular the 
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applications we seek for are to maps from Lorentzian to Riemannian 
manifolds. A C°° map $ : C(M) — > N into a real //-dimensional 
Riemannian manifold (N, h) is biharmonic if $ is a critical point of the 
functional 



z JC{M) 

The tension field r($) is the C°° cross-section in the pullback bundle 
Q^TN C(M) locally given by 



where □ is the Laplace-Beltrami operator of Fg (the wave operator as 
Fg is Lorentzian) and = y % o $. Also if (U, x A ) is a local coordinate 
system on M such that C 7r _1 (C7) and 7 : ^{U) ->• R is a 

local fibre coordinate on C(M) then (7r _1 ([7), w A = x A o 7r, u 2n+2 = 
7) are the naturally induced local coordinates on C(M) and [F pq ] = 
\Fpq\~ 1 while F pq = Fg(d p ,d q ). Here d p is short for d/du v . Finally 
Xf is the natural lift of <9j = d/dy l i.e. the local smooth section in 
$ _1 TiV -> C(M) given by Xf (z) = (3i)* w for any z G S^OO- 
The bundle metric h appearing in is naturally induced by h in 
Q^TN ->■ C(M) so that /i $ (Xf ,X?) = ^ o $. Our main result is 

Theorem 1. Let M be a compact strictly pseudoconvex CR mani- 
fold, of CR dimension n, and 9 a contact form on M with Gg posi- 
tive definite. Let $ : C(M) — )■ N be a smooth S l -invariant map and 
: M — > N the corresponding base map. Then E 2 ($) = 2irE 2 fi{4>)- 
Consequently if $ is biharmonic then is a critical point of E 2 ■ The 
Euler- Lagrange equations of the variational principle S E 2 ^((p) = are 



of N. Consequently the vertical lift to C(M) of any C°° solution f to 
(J3J) is a biharmonic map (with respect to the Fefferman metric Fg). 

Here 71,(0) = r(0; 6, V h ) is the subelliptic tension field of : M — > N 
i.e. the C°° section in _1 TX — > M locally given by 



where A;, is the sublaplacian of (M, 6) (cf. the Introduction or Defi- 
nition 2.1 in [17], p. 134), 0* = y 1 o and Xf(x) = (c^)^ for any 
x G 0~ 1 (V) (the natural lift of d/dy l as a section in <p~ l TN — > M). 



(3) 






r b (0) = { A 6 l + 2g a ? (V jk o 0) T a (0^')T ? (0 fc )} X\ 
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Note that Xf = Xf o n. The rough sublaplacian is the second order 
differential operator locally given by 



2n 



0=1 

for any V G C oo (0 _1 TA r ), where {X a : 1 < a < 2n} is a local orthonor- 
mal (i.e. G e (X a ,X b ) = 5 ab ) frame of H(M). Moreover 0~ 1 V fc is the 
connection in _1 TA r — > M defined as the pullback of the Levi-Civita 
connection V h of (N, h) by <p i.e. 

(r^h/e^xt = ^ (T* k o 0) x*. 

Given a bilinear form B on T(M) we denote by Il^I? the restriction 
of B to # (M) ® H(M). We shall need the following 

Theorem 2. For any smooth map : M — )• iV the rough sublapla- 
cian Af is a formally self adjoint (with respect to the L 2 inner product 
(V,W) = f M h+(V,W)0 A {d0) n , V,W G C^^TN)) second order 
differential operator locally expressed as 

( 2n 

(5) Af V = \ A b V' + 2 W) (T* jk o 0) X a (V k )+ 



0=1 

2n 



(r* fc o 0) a^' + g X a {<t?)X a {f l ) [ + nfi. 



V k V XT 



where V = V i Xf , for any local orthonormal frame {X a : 1 < a < 
2n} of H(M) defined on the open set U C M and any local coor- 
dinate system (V,y l ) on N such that 0~ 1 (V) C U. Let D* be the 
formal adjoint of D = (0~ 1 V /l ) H i.e. (D*ip,V) = (ip,DV) for any 
p G C°°(H(M)* <g> (p^TN) and any V G C%° (^T N) . Then 

(6) Af = —D*D. 

In particular (A b V, V) < 0. 

The proof of Theorem [2] will be given in § 3. If N — M. u then (by (jlj)- 
(|5J)) the subelliptic biharmonic map equations become L<p l = where 
L = A b o A b (the bi- sublaplacian) is a fourth order hypoelliptic operator. 
The analysis of the scalar case N = R (maximum principles, existence 
of Green functions for L, Harnack inequalities for positive solutions to 
Lu = 0, etc.) is however open. The calculation of the Green function 
for A 2 = A o A (where A is the ordinary Laplacian on R") is due to T. 
Boggio, [10]. The existence of the Green function for A b follows from 
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work by J.M. Bony, [12], while estimates (on the Green function and 
its derivatives) were got by A. Sanchez-Calle, [29J, yet the problem of 
adapting their techniques to the bi-sublaplacian is unsolved. 

3. The rough sublaplacian 

The second order differential operator is similar to the rough 
Laplacian on vector fields due to G. Wiegmink, [3D], and CM. Wood, 
[3T] . Let <p : M — > N be a smooth map and <j)~ l TN — > M the pullback 
bundle. Let h? (respectively _1 V h ) be the pullback of the Riemannian 
metric h (respectively of the Levi-Civita connection V h ) by 0. Then 
is parallel with respect to 0~ 1 V /l . We shall establish 

Lemma 1. For any V, W G C oo (0~ 1 TA^) there is a smooth tangent 
vector field X$ on M such that 

(7) h+(At V, W) = A b [h\V, W)\ + h+{V, AfW)-2 div(X^) 

where the divergence is taken with respect to the volume form \I/ = 
9 A (d9) n i.e. C x ^> = div(X^)^. 

Proof. Let {X a : 1 < a < 2n} be a local orthonormal frame of 
H(M). As (0- 1 V h )^ = 

hf>((<f ) - 1 v h ) Xa (<t>- 1 v h ) Xa v J w) = 

= x a (h+ ((fv^/, w)) - h* {{r^ h )x a v , {F 1 v h ) Xa w) = 
= x 2 a {h*(v,w))-2x a {h* (v, (^v h ) Xa w))+h^ (v, (r'v^w), 

h*{(r 1 V h h XaXa V,W) = 
= (\7 Xa X a )(hHV,W)) - h+(V, (r^ h h Xa x a W). 

Let us recall that A b u = div(V H u), u G C 2 (M), where V H u G 
C°°(H(M)) (the horizontal gradient of u) is given by V H u = n#Vu 
and g e (Vu,X) = X{u) for any X G 3C(M). Also U H : T(M) -> H{M) 
is the natural projection associated to the direct sum decomposition 
T(M) = H(M)®RT while g e is the Webster metric of (M, 0) (cf. Def- 
inition 1.10 in [T7], p. 9). A large amount of the existing subelliptic 
theory is built on the Heisenberg group (cf. [TT], p. 155) i.e. on the 
noncommutative Lie group H n e C x 1 with the multiplication law 

(z, t) ■ (w, s) = (z + w, t + s + 2 lm(z ■ w)) 

for any (z,t), (w, s) G H n . Let Z a be the Lewy operators i.e. 

d . d 
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Then Ti t0 (W n ) = Yla=i is a CR structure on H n , of CR dimension 
n, making H n into a CR manifold (cf. e.g. [17J). The Levi distribution 
H(M n ) is spanned by the left invariant vector fields X a G X(H n ) given 
by 

d d d a d 

Xa = d^ + 2y dt' Xn+a = " 2x w 

and the horizontal gradient is familiar (cf. [TT], p. 68) in subelliptic 
theory as the horizontal iJ-gradient i.e. V H u = J2 2 a n =1 X a (u)X a . Note 
also that for the Heisenberg group H = — YldLi -^a (Hormander's sum 
of squares of vector fields) as X* = —X a . 

Let V be the Tanaka- Webster connection of (M,6). As = 
the divergence of a vector field may also be computed as the trace of 
its covariant derivative with respect to V. Hence A b u may be locally 
written as 

2n 

A b u = Y,{Ku-iy Xa X a )u). 

a=l 

Consequently 

(8) h* (A* V, W) = A b [h+(V, W)] + (V, Af W\ + 

a 

Let X^ G H(M) be the vector field determined by 
GeiXfrY) = h+iV, {<f>- x V h ) Y W) 
for any Y G H(M). Then (by Vg e = 0) 

a a 

= y^{ge{V x X f , X a ) + ^ , V ' x X a )} = 

a 

= div(x ) + £V(v, {r x v h )v XaXa w). 

a 

Indeed (^(VtA^, T) = (as H(M) is parallel with respect to V). 
Together with (jSj) this leads to ©. Lemma Q] is proved. 

Let us assume that either M is compact or at least one of the sec- 
tions V, W has compact support. At this point we may integrate (j7j) 
over M and use Green's lemma to show that (Af V, W) = (V, Af W). 
Moreover, if V = V { Xf is a C°° section in (j)~ l TN -»■ M then 

(9) (^ h )xV = + X^)V k (T* h o 0) }X? 
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for any X G X(M). The proof of (jSJ) follows from 

{r^fxV = {x 2 (v l ) + 2xtf) (r jk o /) x(v k )+ 



+ 



v K > XT. 



Let D = (4>- 1 V h ) H i.e. G C°°(H(M)* ® <p- l TN) is the restriction 
of (0" 1 V /l )K to H(M). An L 2 inner product on C°° (H (M)* ® (j)~ l T N) 
is given by 

« In 
J M a=l 

for any 99, ip G C°°(H(M)* (g) <fi~ l TN) and any local orthonormal frame 
{X a : 1 < a < 2n} of #(M) on U C M. Then for any V G 

(d*<p,v) = I J2 h H<px a , [<i>- l v h ) Xa v)t> = 

J a 

= I ^{X a {h*(< P X tt ,V))-h*{(<l>- 1 V h ) Xa <pX a , V)}-*. 

Let X v y G H(M) be determined by G e (X v y, Y) = h^{ V Y, V) for any 
Y G H(M). Then 

a a 

= 52{9e(Vx a X v , v ,X a ) +g e (X VtV , V Xa X a )} = 

a 

= div(X^v) + h+(<pV Xa X a , V). 

a 

We may conclude that 

2n 



a=l 



on [/ and then D*DV = -AfV for any V G C^^TN). 
Proposition 1. The symbol of the rough sublaplacian is 
(10) a 2 [xt)j= [u(T x ) 2 - \\u\\ 2 ] v 

for any u G T*{M) \ {0}, v G {(j)~ 1 TN) x and x G M. Therefore A* ^s 
a degenerate elliptic operator and its ellipticity degenerates precisely in 
the cotangent directions spanned by 9. 
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Proof. Let T'(M) = T*{M) \ (0) and let II : T'(M) — >■ M be the 
projection. If E — > M and F — > M are vector bundles we set 

SmblfcCE, F) = {ae Rom(U*E , U*F) :o pw = p k o^, p > 0} 

(with k G Z). Let <7fc(L) G Smblfc(_E, F) be the symbol of the k- 
th order differential operator L G DM)-(E, F). We wish to compute 
<r 2 (Af) G Smbl 2 (0- 1 TA^, (fi^TN). To this end let w G T'(M) such 
that n(w) = x and let / G C°°(M) such that = uj. Also let 

u G {<p- l TN) x and 1/ G C°°{(f)~ l TN) such that 14 = Then 

* 2 (Af) = ~A{ (*). 

Then (fTOl) follows from the identities 

A 6 (m 2 ) = 2m A 6 m + 2|| V H u\\ 2 , 

AffoV) = g A{V + (A b g)V + 2(r 1 V%» g V , 
where g = u 2 and u — f — f(x). The norm in fTTOT) is ||a;|| = gg x (uj,uj). 



Hence for each v G Ker 
A G E\ {0}. Q.e.d. 



a 2 (A 



6 



either t> = or u> = X8 T for some 



4. The first variation formula 

Let M = M x (-e, e), e > 0, and let F : M ->■ X be a smooth 1- 
parameter variation of by smooth maps i.e. 0o = where <fit = F octt 
and a 4 : M — )■ M is the injection = (x,t) for any x G M. Let 

V G C°°(<j)~ 1 TN) be the corresponding infinitesimal variation 



<9 

V x = (d( x , )F) — 



, x G M. 

(*,0) 

Let us consider the second fundamental form of <fi (cf. R. Petit, [28J 

M)(x, y) = {r 1 v h ) x <p if Y - 0*v x y, x, y g x(m). 

Here 0*X denotes the cross-section in <j)~ 1 TN — )■ M given by 

(<M0(x) = {dx<f>)X x , xGM, 
for each X G X(M). Then (cf. e.g. [8]) 

2n 

r 6 (0) = trace Ge n H /3 6 (0) = ^A,(0)(X o ,X o ) 

a=l 

on [/ C M. We shall establish the following 
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Theorem 3. Let M be a compact strictly pseudoconvex CR manifold 
and (J) : M — >■ N a smooth map into a Riemannian manifold. Then for 
any smooth 1-parameter variation F : M — >■ N of <p 

(11) j t {E2, b (<Pt)} t=0 = J m h* (V, M 6 (0)) vi; 

where R h = </> _1 i? /l is given (jl h (u, v)wj = R^ x ^(u x ,v x )w x , for 
any u,v,w G C oo (0 _1 TA^) and any x £ N. 

If M is not compact one may, as usual, integrate over an arbitrary re- 
latively compact domain Q C M and consider only smooth 1-parameter 
variations supported in Q. 

Proof of Theorem [3J Given X e X(M) we define X E X(M) by 
setting 

X(x,t) — (d x a t )X x , (x,t)eM. 

Then 

(((f>t)*X) x = (d x 4>t)X x = {d( Xjt )F)(d x a t )X x = (d( x , t )F)X( X j) 

hence 

(12) (<k)*X = (F*X) o a t , \t\ < e, 

with the obvious meaning of F*X as a section in F _1 TA r — > M. Simi- 
larly the identities 

(13) Xp=X[oa u l<i<u, 

relate the local frames {Xf* : 1 < i < v} and {X[ : 1 < i < u} in the 
pullback bundles (p^TN M and F~ l TN M respectively Let 
{X a : 1 < a < 2n} be a local orthonormal frame of H(M). We wish to 
compute 

2n 
o=l 

for any \t\ < e. Note that 

(14) v h ) x ^Y = {x(y^) + x(^)y (<p k ) (r;. fc o 0) }xf 

where (jf = y % o <p. As (f)\ = F l o a t it follows that 

(15) Xiti) = X(F*) o a t , X(Yti) = [X(YF*)}oa t , 
for any X, Y G X(M). Therefore (by flUjl-fllS])) 



(16) {<j>^V h ) x {<t>t),Y = {F-^jtFjT 



l-lT-?h\ 



a t 
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where F 1 V h is the connection in F 1 TN — > M induced by V h . Then 
(by (HID and (USD) 

(17) n{<Pt) = n(F) o a t , \t\ < e, 

where r b (F) G C°°(F~ l (V), F^TN) is defined by 

2n 



a=l 



Let /i F be the metric induced by h in F l TN ->Mso that 
h+ (X? ,Xp) = h F {X[ ,Xf)oa t , l<i,j< v. 
Consequently 

In 

\\n(4>t)\\ 2 = J2 hF ({F~ l V h )xF*X a -F*V^X a ,T{F) 



o a t 



a=l 



so that 



-{E2, b (<Pt)} t=0 = l 



dt 



M 



d_ 

di 



{h F (r(F),r(F))} [xfi) *(*) 



2/i 



a=l 



J / (x,0) 

Let i? F lvh be the curvature tensor field of F~ 1 V h . As 



X £ 
a ' dt 



0, 1 < a < 2n, 



it follows that 

(is) (rV")^^-^)^ = 



On the other hand 



(F-V^/afFj - (i^V^F*^ 



<9 



<9 
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so that 
(19) 

By (Jlgjl-flig 



(20) |{^, 6 (^)} t=0 = £ / A / F ((^ 1 V' l )x a (^- 1 V h )i 



(*,0) 



We compute separately the integrand in the right hand side of (!20|) as 
follows. First (as {F~ l V h )h F = 0) 

h F ((F- 1 V h ) jfa (F- 1 V fc ) A .F.| J 7 4 (F)) = 



(( F-1V ")x a ^|' ( J? " 1V *)*. 7 *( F )) 
= A > a X (l (V (F^,r 6 (F) 



5 



-2X a ^ (F*- , (F^V^^F) j j + 

Moreover, as X((p) o a t = X((p o a t ) for any X e Af(M) and any 
V? G C°°(M) 



(21) 



X X a U*'(F*|-,T(F) 



at 



o a t 







X a X a [h" lF*—,n(F))oa t 



dt 



Another tautology we shall need is 



V 



d_ 

"di 



o a . 



The following calculation 

[(^^(f)]^ 

= {X(r b (FY) + X{Fi)r b {F) k (T% o F)} (x>t) Xf (s,t) 
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= {X{n{<p t y) + X{<pl) n {<p t ) k (T) k o <p t )} x xt\x) 

shows that 

(22) [{F~ l V h )^n{F)\ o a t = {<j>T l V h ) Xa n{<t>t) ■ 

Let X t e H(M) be the horizontal tangent vector field determined by 

G e (X t ,Y) = h F ^^(F-^yniF^J oa t 

for any Y G H(M). Then (by = 0) 
d 



Xa { h " ^" 1V/1 )^ T ^ F )JJ oa t = X a {G e (X t ,X a )) = 

= g e (V Xa X t ,X a )+g e (X t V Xa X a ) 

that is 

(23) X a (h F fF,|,(F-'V') 1 T l (F)jj oa t = 9 ,(V x „X t , A'„)+ 



A calculation similar to that leading to (1221) furnishes 

(24) [(F-V^(F-V*) ? t 6 (F)] oa f = (0^)^(0*). 
Summing up the information in (12T1) and (I23l) -( j24l) 

(25) fV ((F-W h ) x SF-W h )x^ t ,r b (F)\oa = 

a=l ^ ' 

= h* (y, A^r fe (0)) -2div(X ) + 
+ £{X fl X fl r 6 (0))) - fc* (V, (0- l V*) Vx .A.74(0))}. 

a 

Using (125|) and 

^((F-V')^F,|,r t (F) 
we may conclude that 

(26) E^l^" 1 ^.^" 1 ^)*-^^- 

a=l ^ 

-(^ lv %^-|' T ' (F) ) oc, " = 
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= h* (y, At r 6 (0)) - 2div(X ) + A b [h+(V, r h (<f>))] . 

Given local coordinates (U, x A ) on M we set x A = x A o p where p : 
M — > M is the natural projection. To compute the curvature term in 
the right hand side of (120]) we conduct 



R F -^ h (X,^Y 



= X A (x)F(0^) x {(F- 1 V ?i ) 9/fe 4^- 1 V /l )^Xf- 
-(F _1 V /l )a/9t(i 71_1 V ?i ) 9 / 9:c AX/ , }( :r> o) = 

: X- 4 (x)F(0^) x .{(F- 1 V ?i ) a/ ^[^ r (1% o F) Xf ] 
-(^V*),,*^ {T%oF)X[]} m = 



x A ( x )Y{<p) a 



9 2 F k (r oF) + ^l ( d Fk oF) ^l\ X F 

dtdx A { kj ' dx A 1 dy e ' dt j 1 



r)F k r)F l 
— ( r ti oF)—(T 7l oF)X^ 



(x,0) 



apt r)W^ 
= X A (x)Y( ( J ) % — (x J 0)^ I (x,0)x 

= X A {x)Y^) x (R% k .{<l>{x)) — {x,U)^{x,U)Xi{x) 
where X = X A d/dx A on U. Noticing that 

X(0%=X A (x)^(x,O), xeU, 



we obtain 

dt 



(27) (i?^ lvh (X, %-)F m Y 



(x,0) 
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= X^VW)*— (x,0) {R h )l, (<P(x))Xt(x) 



Vx =—{x,Q)Xt{x), XeU, 



Finally 
yields 

(28) (V- lv,i (x, ^ = ^((0^0., V X )(<PX) X 

for any a; G M. Next (by ( BHD anc ^ the symmetries of the Riemann- 
Christoffel 4-tensor of (N, h)) 

h F U F - Xvh {X a , ^)F*X a , n{F)) = 

= K(x) (Rl (x) (((j)*X a ) x , V x )((/>JCa)x , n{(f>)x) = 

= -h^x) (R^ x ){n{<j))x , {4>*x a ) x ){<f)*x a ) x , v x ) . 

Together with f )26|) and Green's lemma this leads to the first variation 
formula (ITT)) for i?2,6 in Theorem |3j 

The following concept is central to this paper. A smooth map <fi : 
M — > N is said to be a subeliptic biharmonic map if is a critical point 
of the functional E 2t b : C°°(M, N) — > R. By Theorem [3] a smooth map 
<p : M — > N is subelliptic biharmonic if and only if is a solution to 
(jlj). A pseudoharmonic map is trivially subelliptic biharmonic. 

5. Subelliptic biharmonic maps and Fefferman's metric 

We start by proving the identity E2(0 o tt) = 2ttE2^{4>) in Theorem 
[lj Let us set 3a = Af T B with A^ e C°°(U). As to the range of indices, 
we adopt the convention 

A,B,C,--- e {0,1,- •• ,n,T,-- - ,n} 

with T = T. Then (by ([I])) 

Fab = G e (d A , d B ) + 6(d A )a(d B ) + 6(d B )a(d A ) = 

= QcdWl + ^l) + AVb + \>A 
where a a = 0"(<9a)- A calculation based on (j2J) shows that 

u a = — *— {*Af (T" - -g^T B {g a -B)) - , 9 r A ,} o tt 
n + 2 l AK Ba 2 y 4(n + l) A 
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where T Ba are given by Vr B T a = T Ba Tp. Moreover (by (CD)) 
F A , 2n+2 = 2[(tt*0) a](d A , d/dj) = — |— X° A , 

^2n+2,2n+2 = 0. 

Next, using F ab F bc = 5 a c (with a, 6, c, • • • <E {1, • • • , 2n + 2}) we find 

p AS p . A C wl,2n+2 _ 



(29) 



F AB X° B = 

x° 

p2n+2,Bp _|_ A C ^2n+2,2n+2 _ g 
71 + 2 



F 2n+2,B A j = „ + 2 . 

Taking into account that d& /d'y = and F AB X° B = we have 
(30) jm^) 



F AB (T% o $) j AIA^T Q (^)T^(0 fc ) + AlAjT a (^)T-^^ 



+ X a A X B T^)T^ k ) + A^W')T^ fc ; 



O 7T. 



We need the following 

Lemma 2. T/ie (reciprocal) Fefferman metric is related to the (recip- 
rocal) Levi form by 

(31) F AB X a A x\ = tf*, 

(32) F AB X a A X B = 0. 

Proof. The identities (129]) may be written 
^ AS ^(ASAg + AgA|) + F ab X° c (Jb + ^ i^-+ 2 A c = 

F AB X° B = 0, 

, 77>2n+2,B \ C> , ^ 772n+2,2n+2 \ n 

+^ A c ctb + ~ - 2 J = U, 

F 2n+2 ' B X° B = n + 2. 
If := A -1 then (by the first of the previous four identities) 

A / I r-iA.2n+2 | t?AB _ \c0 , 77MB /\«rf+n , \^ra\ 

= CrXo + aB >° D + 9ap{ X B° D + X B b D ) 

ill \~ £ 
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yielding 



(33) 



( 




pA,2n+2 _|_ pAB 



The second and third of the identities (1331 lead to ( 13TT) and ( 1321) . re- 
spectively. Lemma |2] is proved. 

By Lemma [2] we may write (1301) as 



(34) F pq (r jk o $) — — = 2 {(r;. fc 0) 



By a result of J.M. Lee (cf. [14], or Proposition 2.8 in [T7], p. 140) 

□ ( U O tt) = (A 6 U) O 7T, iiGC 2 (M). 

Hence (by (|34|)) 



On the other hand if 7r _1 ({7) ^ [/ x 5 1 is a local trivialization chart of 
the canonical circle bundle and u G C°°(M) is a function supported in 
U then 



(by integration along the fibres of S 1 — > C(M) — > M, cf. e.g. (2.49) in 
[IT] , p. 141). Hence (by a partition of unity argument) 



To prove the next statement in Theorem [T] let {4>t}\t\<e be a smooth 
1-parameter variation of (0o = 0) so that $ t = t o7r is a 1-parameter 
variation of $. Therefore (by (|36l) ) if $ is biharmonic then is a critical 
point of E 2 fi- 

The converse doesn't follow from fl36l) but rather from the first vari- 
ation formula for i?2,&. Indeed let : M — > N be a smooth solution 
to (jl]). A slight modification of Jiang Guoying's arguments (cf. [2"2"] ) 
leads to the following 

Lemma 3. A smooth map $ : C(M) —tNis biharmonic if and only 
if $ is a solution to 



(37) = □ $ r($) + tracer {($- 1 i? /l )(r($), $* •)$„■} = 



(35) 





(36) 
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where 

2n+2 

(38) □*» = J2 <*{(* -1 V% - (^V*)^, }9J 

p=i 

the rough Laplacian on C(M). Here 53 G C°°($ _1 TA^) and {X p : 
1 < p < 2n + 2} is a local Fg-orthonormal (i.e. Fg(X p ,X q ) = e p 5 pq with 
e l = ■■■ = e 2n+ i = -e 2n+2 = 1) frame in T(C(M)). Also V C(M) is the 
Levi-Civita connection of (C(M), Fg). 

Let {X a : 1 < a < 2n} be a local Gg-orthonormal frame in H(M), 
defined on the open set U C M. Then {Xj, ± S : 1 < a < 2n} is 
a local Fg-orthonormal frame of T(C(M)). Here for any X G X(M) 
we denote by G X(C(M)) the horizontal lift of X with respect 
to the connection 1-form a on the canonical circle bundle (thought of 
as a principal ^-bundle over M). We recall that X\ G Ker(<r) z and 
(4vr)X^ = X 7r(z) for any z G C(M). Also 

n + 2 9 

b ~ 2 9 7 • 

Let 1/ G C oo (0- 1 TJV) and = V o 7T. Let $ = o tt. The identities 

(i-^wtf = o, 

(<*> - l V h ) d/duA Xt = [(r^a/a^Xt] o tt, 

imply 

(39) ($ _1 V V93 = K^'V^x^] o tt, 

(40) ($- 1 V ?i ) x t($~ 1 V /l ) y t2J = [(^ _1 V h )x(^" 1 V' l )yV] o TT, 
for any X, Y G T(M). At this point we need 

Lemma 4. For any X, F G H(M) 

V C X \ M) Y^ = (VxY? - (d0)(X, F)T^ - L4(X, Y) + (d<r)(X^ F^S, 

V ; ( t M) Tt = (rX + ? X)t, 

V^ M) X^ = (V T X + gX) f + 2(da)(X\T t )S, 

V C X [ M) S = V? M) Xt = (JX)\ 

V C T \ M) T^ = Q\ vf M) S = ^ 
v c ( M) Tt = v c(M) 5 = 0; 

w/iere q : H(M) ->■ H(M) is given by G 6 (qX,Y) = (da)(X\Y r ) and 
Q G H(M) is given by Gg(Q,Y) = 2(da)(T t ,Y r ) . Also r t/ie psen- 
dohermitian torsion o/V and A(X, F) = Gg(rX, Y). 
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Cf. Lemma 2 in E. Baletta et al., [9], p. 083504-26. As a consequence 
of Lemma H] 

Vj ( t M) X t = (V x X) t - A(X,X)S, 
V^I* ± 5) = Q\ 

hence (by (I39|)) 

a x <> 

= E^^W^t - A (*a, X a )($- 1 V ft ) 5 }QJ = 
a 

= J2l^ h h Xa x a V]on 

a 

as trace^A = trace(r) = (cf. e.g. (1.59) in [T7], p. 37). Together 
with flU and 

($- 1 V h ) s 9J = 
this allows one to conduct the following calculation 

2n 

□*QJ = J]{($- 1 V' 1 )^ - ($- 1 V /l ) v c ( M, x; r}QJ+ 
+ {($- 1 V A )^ t+5 - ($- 1 V ft ) v a (M ) Tt+s }QJ- 

V T T+S 

rt-s 

a 

that is 

(41) = ( Af V) O 7T. 

It remains that we compute the curvature term in (1371) . As $*X^ = 
(4>*X) o vr for any X, Y e *(M) 

($- 1 i? /l )(r($),<l>,A t )<l>,r t = [(0- 1 i2 fc )(7i(^) J 0,X)^y] OTT 

hence 

tracer {($- 1 i? /l )(r($), } = 

= E($- 1 ^)(r($),$,X:)$,Xt+ 

a 

+ ($- 1 i? h )(r($), ^(T 1 " + S))$.(T* + S)- 
-($- 1 i? /i )(r($), _ S , ))$*(T t - 5) = 

= trace G() 7r H {(^ 1 i? ?i )(r 6 (0), 0* }+ 
+2{($- 1 J R /l )(r($), + ($- 1 J R /l )(r($), $^)$,T^} 
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and $*S f = so that 

(42) tracer {($- 1 i? /l )(r($), ■)$* ■} = 

= trace G() vr^{(0~ 1 J R ?i )(r b (0), 0, •}■ 
Finally (by ()4ip - (|42p ) <p o it is biharmonic. Theorem [1] is proved. 

6. Conclusions and open problems 

We introduced the new concept of a subelliptic biharmonic map as 
a C°° solution to the system (j3J). This is a quasilinear system of vari- 
ational origin whose principal part is the bi-sublaplacian A 2 . A 2 is a 
fourth order hypoelliptic operator, though not elliptic, so that our work 
is part of the program outlined in [23] . Higher order degenerate elliptic 
equations, say of the form A^u = 0, were not studied in the present 
day PDEs literature (cf. e.g. jl] for the elliptic case). Our main re- 
sult is a geometric interpretation of subelliptic biharmonic maps within 
Lorentzian geometry i.e. each C°° solution to (jl]) may be characterized 
as the base map corresponding to a ^-invariant biharmonic map from 
the total space C(M) of the canonical circle bundle endowed with the 
Fefferman metric. Although, as shown in § 4, it makes sense to look for 
weak solutions to (J3J) our methods in this paper are purely geometric 
and a study of local properties of weak subelliptic biharmonic maps ap- 
pears nowhere in the mathematical literature. Neither may the partial 
regularity theory be naively reduced to that of ^-invariant biharmonic 
maps, as C(M) is Lorentzian so that no natural distance function on 
C(M) is available a priori. If Q C M is a bounded domain the func- 
tional E 2t b((j>) and f n |A fe 0| 2 \1> (as introduced by S-Y.A. Chang & L. 
Wang & P.C. Yang, [T5], in the elliptic case for maps : — >• S u ) have 
not been compared so far (and of course the corresponding regularity 
for S J n |A;,0| 2 \1/ = is unknown). 
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